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Abstract  In  the  fast  time-varying  shallow  water  environment,  the  phone-space 
adaptive  processing  using  sample-matrix-inversion  (SMI)  approach  outperforms 
the  iterative  least-mean-squared  (LMS)  approach  due  to  its  rapid  convergence. 
The  SMI  approach  uses  singular-value-decomposition  (SVD)  to  decompose  a 
sample-matrix  into  a  set  of  spatial  eigenvectors  and  their  associated  eigenvalues. 
Adaptive  beamforming  is  applied  by  matching  the  steering  vectors  with  the 
spatial  eigenvectors  weighted  by  their  eigenvalues.  When  a  long  towed-array 
undergoes  significant  maneuvering,  the  array  shape  and/or  the  target-to-array 
relative  geometry  changes  rapidly  within  a  processing  interval.  A  signal  in  the 
phone-space  sample-matrix  is  split  into  more  than  one  eigenvector  resulting  in 
signal  mismatch  in  the  subsequent  beamforming.  Its  power  is  split  into  many 
eigenvalues  resulting  in  signal  loss  in  the  processing. 

In  subarray  beam-space  adaptive  beamforming  (SABS  ABF),  a  beam-space 
sample-matrix  is  formed  at  each  search  cell  by  focusing  subarrays  to  a  cell 
using  the  dynamically  updated  array  shape  at  each  time  step.  The  SMI  adaptive 
beamforming  then  is  done  by  decomposing  the  beam-space  sample-matrix  and 
matching  the  beam-eigenvectors  with  a  unity  steering  vector.  The  beam-space 
sample-matrix  has  a  lower  rank  than  the  phone-space  sample-matrix  so  that  a 
stable  estimation  can  be  reached  with  fewer  time  samples.  The  dynamic  array 
shape  compensation  in  SABS  ABF  makes  a  signal  in  the  beam-space  sample- 
matrix  less  likely  to  be  split  by  the  SVD.  The  subsequent  processing,  after  the 
SVD,  experiences  less  signal  mismatch  and  signal  loss.  A  dynamic  cable  model 
is  used  to  simulate  a  long  towed-array  going  through  turns  in  a  measured  shallow 
water  current  field.  Simulations  show  that  significant  signal  loss  in  phone-space 
SMI  processing  is  recovered  in  SABS  ABF  processing. 
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Abstract: 

In  the  fast  time-varying  shallow  water  environment, 
the  phone-space  adaptive  processing  using  sample- 
matrix -m  vers  ion  (SMI)  approach  outperforms  iterative 
least-mean- squared  (LMS)  approach  due  to  its  rapid 
convergence.  The  SMI  approach  uses  singular-' value- 
decomposition  (SVD)  to  decompose  a  block  sample- 
matrix  into  a  set  of  spatial  singular  vectors  and  their 
associated  singular  values.  Adaptive  beamforming  is 
applied  by  matching  the  steering  vectors  with  the  spatial 
singular  vectors  weighted  by  their  singular  values.  When 
a  long  to  wed-array  undergoes  significant  maneuvering,  its 
shape  relative  to  a  target  changes  rapidly  within  a 
processing  interval.  The  target  signal  in  the  phone-space 
sample -matrix  is  split  into  more  than  one  singular  vector 
resulting  signal  mismatch  in  the  subsequent  beamforming. 
Its  power  is  split  into  many  singular  values  resulting 
signal  loss  in  the  processing. 

In  subarray  beam- space  adaptive  beamforming 
(SABS_ABF),  a  beam-space  sample-matrix  is  formed  at 
each  search  cell  by  focusing  subarrays  to  a  cell  using  the 
dynamically  updated  array  shape  at  each  time  step.  The 
SMI  adaptive  beamforming  then  is  done  by  decomposing 
the  beam- space  sample-matrix  and  matching  the  beam- 
eigenvectors  with  a  unity  steering  vector.  The  beam- 
space  sample -matrix  has  a  lower  rank  than  the  phone- 
space  sample-matrix  so  that  a  stable  estimation  can  be 
reached  with  fewer  tune  samples.  The  dynamic  array 
shape  compensation  in  SABS  ABF  makes  a  signal  in  the 
beam- space  sample -matrix  less  likely  to  be  split  by  the 
SVD.  The  subsequent  processing,  after  the  SVD, 
experiences  less  signal  mismatch  and  signal  loss.  A 
hydrodynamic  cable  model  is  used  to  simulate  a  long 
towed -array  going  through  turns  in  a  measured  shallow 
water  current  field.  Simulations  show  that  significant 
signal  loss  in  the  phone- space  adaptive  processing  is 
recovered  in  the  SABS  ABF  processing. 

1.  Introduction 

There  are  three  mam  issues  for  a  long  towed-array 
going  through  dynamic  maneuvering,  they  are  increasing 
of  array  self-noise,  rapid  changing  in  array  shape  (array 
element  location  AEL),  and  signal  spreading  in  the 
phone-space  sample -matrix.  Premus  et.  al.  [1]  showed 


the  use  of  a  beam-space  adaptive  beamforming  for  towed 
array  self-noise  cancellation.  Ianniello  et.  al.  [2] 
demonstrated  array  shape  estimation  through  turns  using  a 
hydrodynamic  model  with  data  from  the  non-acoustic 
heading  and  depth  sensors.  This  paper  addresses  the  issue 
of  signal  spreading  that  the  received  signal  on  an  array 
changes  within  the  processing  interval  due  to  array 
motion. 

Figure  1  shows  the  simulation  scenario  using  in  this 
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Figure  1:  Simulation  Geometry 

study.  A  towed  ship  heads  toward  North  at  the  beginning 
of  the  simulation,  it  makes  an  180°  turn  in  10  minutes 
then  heads  toward  South.  A  long  array  with  210 
hydrophones  and  a  spacing  of  3.14  m  is  towed  behind  the 
ship  with  the  first  phone  at  360  m  from  the  ship.  A 
hydrodynamic  cable  model  is  used  to  simulate  a  realistic 
motion  of  the  array  going  through  the  turn  in  a  measured 
shallow  water  current  field.  A  target  starts  at  a  distance  of 
20  km  and  approaches  the  ship  from  East  at  a  speed  of  4 
kts.  In  additional  to  the  own- ship,  two  other  interferers 
are  included  to  simulate  a  stressful  environment.  One 
interferer  starts  from  North-West,  it  crosses  the  back- 
beam  of  the  target  and  heads  to  South.  The  other  starts 
from  South-East  and  moves  into  the  near-field  of  the 
array.  The  simulation  is  done  at  200  Hz  with  a  sampling 
rate  of  every  half  second.  All  processing  in  this  study 
uses  a  time  interval  of  30  seconds.  True  bearing  is  used 
for  displaying  the  beamforming  results. 
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Let  xn(t)  be  the  received  time  series  of  the  nth 
hydrophone,  the  signal  coherence  (SC)  relative  to  the  first 
hydrophone  is  defined  as 


SCR(t) 


K- 1  z 

y  x[{t  +  kAt)  .  xn(t  +  kAt) 

_ k= 0 _ 

+  kAt)\2  +  kAtrf 

k  = 0  £=0 


here  ^=60  and  At=0.5  s.  For  target  only,  Figure  2  shows 


Figure  2:  Time  histories  of  signal  coherence  relative  to 
phone  one 


the  time  history  of  signal  coherence  for  hydrophone  in  the 
middle  of  the  array  (n=105)  and  at  the  end  of  the  array 
( 11=21  (Irrespectively.  The  signal  coherence  is  close  to 
perfect  across  array  before  the  turn  and  is  significantly 
degraded  in  the  turn.  The  effect  of  the  residual  tail 
motion  after  the  turn  is  seen  in  hydrophone  210. 

Let  X  be  a  phone-space  sample-matrix  for  a 
processing  interval  KAt  and 

jq(?)  xl(t-\-At)  xl(t+2M)  ...  .^(£+(£-1)#) 
X^t)  X 2(t+bt)  x2(t-\-2M)  ...  X2(t+(K-l)/st) 


xN(t)  xN(t+tf)  xN(t+ 2M)  ...  XN(t+{K-l)tf)_ 

After  smgular- vahie- decomposition  (SVD),  X  can  be 
written  as 


X=V+AU 


where  V  and  U  are  orthogonal  matrices,  and  A  is  a 
diagonal  matrix 


A  = 


4^ 


0 

0 


°  °  VLT 

The  columns  v,  and  U\  of  V  and  U  are  called  the  left  and 
right  smgular  vectors,  respectively,  and  the  diagonal 
elements  X\  of  A  are  called  the  singular  values. 


Physically,  each  v,«  represents  a  spatial  structure  of  a  signal 
observed  across  the  array  of  hydrophones,  the  associated 
Uj  contains  the  information  of  temporal  variation  of  this 
signal,  and  A;  is  the  mean  power  of  the  observed  signal. 
For  target  only,  Figure  3  shows  the  singular- value  tune 
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Figure  3:  Singular-Value  tune  histories  of  the  phone- 
space  sample-matrix  for  target  only 

histories.  It  shows  that  target  energy  is  contained  within 
one  singular  value  before  the  turn  and  is  spread  into 
several  singular*  values  in  the  turn.  Without  an 
appropriate  signal  model  to  account  for  the  spread,  target 
estimation  using  the  phone-space  sample-matrix  will  be 
significantly  degraded. 

Subarray  beam-space  processing  is  a  two-stage 
processing.  In  the  first  stage,  a  phone-space  sample- 
matrix  is  transformed  into  a  beam-space  sample-matrix 
using  the  dynamically  updated  AEL.  The  transformation 
compensates  the  temporal  variation  of  the  received  signal 
on  the  array  due  to  rapid  changes  in  array  shape.  The 
signal  in  the  beam-space  sample-matrix  is  more  stationary 
and  experiences  less  spread  so  that  the  subsequent 
processing  can  achieve  the  optimal  result.  This  paper  is 
organized  as  follows:  Section  2  reviews  the  white-noise- 
constrained  adaptive  beamforming,  Section  3  describes 
the  subarray  beam- space  adaptive  beamforming 
(SABS  ABF),  Section  4  shows  the  simulation  results,  and 
Section  5  summaries  the  study. 

2.  White-Noise-Constrained  (WNC)  Adaptive 
Beam  forming 

Adaptive  processing  uses  the  measured  signal  plus 
noise  data  vectors  to  minimize  the  sidelobe  contributions 
from  those  components  that  do  not  match  with  the 
steering  vector  for  a  given  search  cell.  Let  R  be  the 
covariance  matrix  of  the  received  signal  and  noise,  and 
R=<XX+>  where  <>  denotes  ensemble  average  over  a 


number  of  sequential  data  vectors  X  and  “+”  denotes 
Hermitian  transpose.  The  minimum  variance 
distortionless  response  (MVDR)  method  minimizes  the 
variance  at  the  output  of  a  linear  weighting  of  the 
hydrophone  array  subject  to  the  distortionless  constraint 
that  signals  in  the  steering  direction  have  unity  gam.  The 
formulation  minimizes  the  variance  given  by 

*W=W+RW 

with  respect  to  the  weighting  W,  subject  to  the  unity  gain 
constraint 


W+A  =  1 

where  A  is  a  steering  vector.  The  MVDR  weight  vector 
W mvdr  can  be  derived  as 


w 


MVDR 


R-1A 

A+R_1A 


where  R'1  is  matrix  inversion.  The  MVDR  output  is 

5  -  1 

mm  A  +  R-iA 

Applying  the  eigen-value  decomposition,  the  covariance 
matrix  R  can  be  decomposed  into  a  set  of  eigenvectors  V,- 
associated  with  a  set  of  eigenvalues  A/,  so  that 

/=l 

and  the  inverse  of  the  covariance  matrix  is  given  by 

N  1 

R_i = y  — v,v,+ 

i= l  A: 

The  MVDR  weight  vector  becomes 

N 


W, 


E|(v;a)v,. 

/=!  2L _ 

Zy  V>+A 

/=l  A: 


The  MVDR  output  becomes 

/  N 


1 MVDR 
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V/=l 


V,+A| 


A"1 


An  alternative  way  to  obtain  V/  and  A/  described  in 
Section  1  is  applying  singular- value  decomposition  to  a 
sample -matrix. 

Without  mismatch  between  data  and  model  and  when 
the  signal  is  loud  enough,  the  steering  vector  A  is 
perfectly  matched  with  signal  eigenvector  in  the  steering 
direction.  The  rest  of  eigenvectors  are  orthogonal  to  the 
steering  vector  and  have  no  effect  on  the  output.  But, 
when  mismatch  is  present,  the  steering  vector  is  no  longer 
orthogonal  to  the  rest  of  eigenvectors.  The  noise  vectors 
associated  with  the  least  significant  eigenvalues  then 
dominate  the  inverse  processing  and  degrade  the  signal 
estimation. 


All  forms  of  mismatch  are  either  deterministic  or 
random.  Deterministic  mismatch  degrades  the  signal 
estimation  and  causes  bias  in  estimation,  but  it  can  be 
minimized  if  a  prior  knowledge  is  included  in  modeling 
the  signal.  Random  mismatch  degrades  the  signal 
estimation  but  will  not  bias  the  estimation.  Random 
mismatch  cannot  be  minimized  so  that  robust  algorithms 
were  developed  to  tolerate  it  to  a  certain  level.  The  white - 
noise-gam  constrained  (WNC)  method  referred  to  by  Cox 
[3],  dynamically  adjusts  the  sensor  noise  level  by  adding 
white  noise  power  to  the  diagonal  elements  of  the 
covariance  matrix  subject  to  an  inequality  constraint  on 
the  sensor  noise  gain.  Adding  noise  expands  the  noise 
space  and  effectively  eliminates  small  eigenvalues  that 
would  otherwise  dominate  the  sum  in  the  MVDR  output 
calculation  due  to  mismatch. 

The  MVDR  white  noise  processing  gam,  defined  as 
the  amplitude  squared  of  the  weight  vector  |WMVDR|2,  is 

directly  proportional  to  the  signal-to-noise  ratio  (SNR). 
The  MVDR  signal  degradation  due  to  mismatch  is 
inversely  proportional  to  the  SNR.  At  low  SNR,  the 
white  noise  processing  gam  approaches  unity  (the 
conventional  linear  processing  noise  gain)  and  the 
mismatch  effect  is  negligible.  At  high  SNR,  the  white 
noise  processing  gain  is  high  and  the  MVDR  output  is 
very  sensitive  to  mismatch.  For  each  steering  direction, 
the  WNC  method  dynamically  adjusts  the  sensor  noise 
level  by  adding  white  noise  to  the  diagonal  elements  of 
the  covariance  matrix.  Adding  white  noise  lowers  the 
“apparent”  SNR  so  that  the  processing  becomes  less 
sensitive  to  mismatch.  Adding  white  noise,  in  the  amount 
8,  to  the  diagonal  elements  of  the  covariance  matrix  is  the 
same  as  adding  8  to  each  eigenvalues  without  modifying 


the  eigenvectors.  The  WNC  weight  vector  Ww/vc  that 
result  is 


w  = 

yyWNC  ft 


E— (v,*a)v„ 

N  1  |2 

f - V,.+A 

tM+s  1  1 


and  the  WNC  output  is 

N  0 

X — '  si: 


trrf. +g)2 


V,+A 


JWNC 


N  1 
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There  are  many  approaches  to  determine  the  amount  of 
white  noise  to  add  in  the  WNC  processing.  In  this  study 
the  following  algorithm  is  used.  For  each  steering 
direction,  the  MVDR  white  noise  processing  gain  and  the 
MVDR  output  are  calculated.  If  the  white  noise 
processing  gain  falls  below  a  pre-selected  constraining 


value,  the  MVDR  processing  is  used.  If  the  white  noise 
processing  gain  is  above  the  constraining  value,  an 
amount  of  white  noise  that  equals  the  MVDR  output  is 
added  in  the  WNC  processing. 

In  this  paper  all  adaptive  beamforming  (ABF)  refer  to 
WNC  adaptive  beamforming  with  a  3-dB  constraint.  The 
phone- space  adaptive  beamforming  (PS_ABF)  refers  to 
the  WNC  adaptive  beamforming  performing  on  the 
hydrophone  sample -matrix  with  a  calculated  phone 
steering  vector. 

4.  Subarray  Beam-space  Adaptive  Beamforming 
(SABS_ABF) 

Subarray  beam- space  adaptive  beamforming  is  a  two- 
stage  processing.  It  has  been  studied  [4]  with  emphasis 
on  its  advantage  of  rank  reduction  for  fast  adaptation,  hi 
this  study,  our  emphasis  is  on  using  the  updated  AEL  to 
compensate  the  temporal  variation  of  a  received  signal 
due  to  rapid  changes  in  array  shape.  Figure  4  shows  the 
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Figure  4:  Subarray  configurations  used  in  this  study 


full- array  and  four  Subarray  configurations  that  are 
considered  in  this  study.  The  notations  in  each 
configuration  are  SEQuential,  number  of  hydrophones  in 
each  subarray,  subarray  overlap,  and  number  of 
subarrays.  For  example,  SEQ_15_10_40  is  a 
configuration  of  15 -phone  subarrays  in  a  sequential 
format  with  overlap  of  10  phones  and  a  total  of  40 
subarrays. 

Let  {7, 2,..., S}  be  a  set  of  hydrophones  in  one  of  the 
subarrays,  the  range-focus  steering  vector  of  the  m* 
subarray  at  search  cell  r  is  calculated  as 


A"<r,')=3sL' 


where 


Vm,!  (0  =  |r  -  rm  s  (r)|  /  c 


c  is  the  sound  speed  and  rms(t)  is  the  measured  position  of 
the  8th  hydrophone  in  the  mth  subarray  at  time  t.  In  the 
first  stage  SABS_ABF  performs  subarray  conventional 
beamforming  using  the  updated  subarray  steering  vectors 
and  transforms  the  phone- space  sample -matrix  into  a 
beam- space  sample -matrix 

bx(j,t)  bfoit  +  At)  ...  bx(j  ,t  +  (K -X)At) 

b2(r,t)  b2(r,t  +  At)  ...  b2(r,t  +  (K-l)At) 

(r  ,0=  ...  . 

_bM(r,t)  bM(r,t+M)  ...  bM(r,t+(K-l)&t)_ 

where 


bjv,t)=  +(r  ,t).  (0 

is  the  conventional  beamformer  response  of  the  mth 
subarray  and  Xm(t)  is  the  data  vector  of  the  mth  subarray  at 
tune  t.  In  the  first  stage  SABSABF  also  transforms  the 
phone-space  steering  vector  into  a  beam-space  steering 
vector.  Let 

A;<r./)=tS>,r)  -  AVr.o] 

where 


=  A^(r.0*Am(r,0 

the  beam-space  steering  vector  is  calculated  as 

Ab(r,t)  =  Ab(r,t)/Ab(r,t) 


For  a  constant  amplitude  non-fading  signal,  the  beam- 
space  steering  vector  is  simply  a  unity  vector  with  a 

uniform  weight  of  1 / V M  on  each  subarray.  In  the 
second  stage  SABSABF  performs  white-noise- 
constrained  adaptive  beamforming  described  in  Section  3 
with  the  calculated  beam-space  sample-matrix  and  the 
beam-space  steering  vector. 


5.  Simulation  Results 


The  simulation  geometry  has  been  described  in 
Section  1  and  a  0-dB  white  noise  is  injected  in  all 
simulations.  For  target-only,  Figure  5  shows  the  WNC 
ABF  responses  along  the  target  track.  Assuming  array  is 
straight  (No  AEL),  PS  ABF  experiences  significant 
mismatch  loss  during  turn.  When  the  measured  mean 
AEL  is  used,  the  mismatch  loss  in  PS  ABF  is  reduced, 
but  there  still  is  more  than  5  dB  signal  spread  loss  in  the 
turn  due  to  rapid  changes  in  array  shape  within  the  30-s 
processing  interval.  SABF_ABF  shows  success  in 
holding  the  target  through  the  turn.  For  target-only, 
Figures  6  and  7show  the  bearing-time  responses  (BTRs) 
resulted  from  PS_ABF  and  SABS_ABF  for  a  0-dB 
target.  The  results  show  that  SABS_ABF  not  only 
reduces  the  signal  spread  loss  but  also  reduces  the 
angular  spread  of  the  target  during  the  turn. 


Figure  5:  WNC  ABF  responses  along  target  track, 
black  curve  is  result  of  PS_ABF  without  AEL,  blue 
curve  is  result  of  PS_ABF  with  mean  AEL,  and  red 
curve  is  result  of  SABS_ABF  for  SEQ_15_10_41 
configuration. 
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Figure  6:  BTR  of  PS_ABF  for  a  0-dB  target 


Figure  7:  BTR  of  SABS_ABF  with  SEQ_15_10_41 
configuration  for  a  0-dB  target 


Figures  8  and  9  show  the  WNC  ABF  results  for  a 
weak  target  in  a  stressful  environment.  In  this  case,  the 
target  power  is  set  at  -18  dB,  the  own-ship  and  the  other- 
two  interferes  are  set  at  15  dB.  The  PS_ABF  result 
shows  that  target  is  lost  in  the  interference  in  turn.  The 
SABS_ABF  result  shows  a  continuous  target-holding 
through  turn  under  such  stressful  condition. 


Figure  8:  BTR  of  PS_ABF  for  a  weak  target  in  a 
stressful  environment 
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Figure  9:  BTR  of  SABS_ABF  with  SEQ_15_10_41 
configuration  for  a  weak  target  in  a  stressful 
environment 


Figures  7  and  9  are  results  of  SABS  ABF  with 
SEQ_15_10_40  subarray  configuration  where  two 
consecutive  subarrays  have  a  66%  overlap  and  the  centers 
are  shifted  by  a  distance  of  5 -phone  spacing.  Figure  10 
shows  result  of  SABS_ABF  with  SEQ_15_0_14 
configuration  where  two  consecutive  subarrays  have  no 
overlap  and  the  centers  are  shifted  by  a  distance  of  15- 
phone  spacing.  It  shows  that  overlapping  the  subarrays 
provides  better  suppression  of  signal  sidelobes.  Figure  11 


shows  result  of  SABS_ABF  with  SEQ_30_10_37 
configuration  where  two  consecutive  su  bar  rays  have  a 
66%  overlap  and  the  centers  are  shifted  by  a  distance  of 
10-phone  spacing.  Signal  sidelobes  are  still  noticeable  in 
this  case.  Figure  12  shows  result  of  S  ABS  ABF  with 
SEQ_30_5_37  configuration  where  two  consecutive 
subarrays  have  a  83%  overlap  and  the  centers  are  shifted 
by  a  distance  of  5 -phone  spacing.  Signal  sidelobes  are 
well  behaved  as  seen  in  SEQ_15_5_40.  This  says  that  the 
distance  between  the  centers  of  two  consecutive  subarrays 
is  an  important  parameter  for  SABS  ABF  processing. 

-13 
-14 
-15 
-1G 
-17 
-18 

-19 
-20 
-21 
-22 
-23 

0  50  100  150  200  250  300  350 

True  Bearing  (deg) 


Figure  10:  BTR  of  SABS_ABF  with  SEQ_15_0_14 
configuration  for  a  weak  target  in  a  stressful 
environment 


Figure  11:  BTR  of  SABS_ABF  with  SEQ_30_10_19 
configuration  for  a  weak  target  in  a  stressful 
environment 
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Figure  12:  BTR  of  SABS_ABF  with  SEQ_30_5_37 
configuration  for  a  weak  target  in  a  stressful 
environment 


6.  Summary 

In  this  paper  we  address  the  issue  of  signal  spread  loss 
in  a  phone-space  sample-matrix  due  to  rapid  changes  of 
array  shape  within  a  processing  interval  when  array  is  in 
maneuvering.  The  loss  significantly  affects  detection  and 
tracking  of  a  weak  target.  SBAS_ABF  uses  the  updated 
AEL  to  compensate  the  temporal  variation  of  the 
received  signal  so  that  signal  in  the  beam-space  sample- 
matrix  becomes  stationary  during  array  maneuvering.  It 
is  shown  that  SABS_ABF  not  only  has  the  advantage  of 
reducing  rank  for  fast  adaptation,  but  also  improves 
target  detection  and  tacking  through  severe  array 
maneuvering. 
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SUMMARY 


Subarray  Beam-Space  ABF  with  updated  AEL  compensation 
recovers  signal  spreading  loss  and  provides  continuous  target 
holding  through  array  maneuvering. 

Beam-space  ABF  involves  extensive  computation 

-  Try  Fast  ABF  such  as  Multi-Stage  Wiener  Filter 

What’s  the  rule-of-thumb  for  subarray  selection? 

-  Distance  between  two  consecutive  subarray  centers  is  an  important 
parameter. 

Is  there  significant  benefit  for  ABF-ABF  vs  CBF-ABF? 


